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TP #1, Ql 

% Nfft: nombre de points de la DFT (fonction fft) 

Nfft = 200000; 

deltaf = 2/Nfft 

xx = [-1 : deltaf : 1-deltaf] ; 

yy = 20*logl0 (abs (f f tshif t (fft (data, Nf f t) ) ) ) ; 
figure , plot(xx,yy); 
figure , plot(xx,yy); 

xlabel ( 1 \omega \ times \pi 1 ) ; y label ( ' dB 1 ) ; 
title (' |X(\omega) | A 2 ' ) ; 

freqz (RIF_1,1) ; 
freqz (RII_b ,RII_a) ; 
figure , zplane (RIF_1 , 1 ) ; 
figure, zplane (RII_b,RII_a) ; 
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TFD du signal audio (module^ 2) 
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Resolution en frequence vs precision 



data = datal + (0 . l*cos (0 . 55*pi* (1 : 40100) ))' ; 

% Le signal a 40100 points; la resolution en 
frequence est determinee par le nombre de points du 
signal et la fenetre qui est utilisee 

Resolution = demi-largeur du lobe principal 

= 2n/40100 = 5xl0 -5 n pour la fenetre rectangulaire 

% A) TFD (fft) calculee sur 200000 points; il y a 
done un point de la DFT a chaque 2n/200000 = 10 -5 n 

% B) TFD (fft) calculee sur 50000 points; il y a 
done un point de la DFT a chaque 2n/50000 = 4xl0~ 5 n 
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L = 40100, fen. Rectang, Nfft = 200000 
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L = 40100, fen. Rectang, Nfft = 200000 
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TP #1, Ql 



figure , zplane (RIF_1 r 1 ) ; 
freqz (RIF_1 ,1) ; 

figure, zplane (RII_b,RII_a) ; 
freqz (RII b,RII a) ; 



Filtre RII 
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Filtre RII ordre 10 
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Filtre RIF 1 




Filtre RIF avec 102 coefficients 
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Filtre RIF 2 




Filtre RIF avec 22 coefficients 
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TP #1, Ql (suite) 

% filtrage du signal musical 
% filtre RIF 

data_firl = conv (data ,RIF_1) ; 
% filtre RII 

data_iir = f ilter (RII_b,RII_a,data) ; 

% realisation du filtre notch 
zl = exp (li*0 . 55*pi) ; 
z2 = conj (zl) ; 
zz = [zl, z2]; 

rr = 10 A ((loglO (0.001) ) /1000) ; 
pp = r r . * z z ; 

[bb,aa] = zp2tf (zz 1 ,pp , 1) ; 
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Poles et zeros 
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Reponse en frequence 





Reponse a I'impulsion du filtre notch (normalisee) 
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Reponse a I'impulsion du filtre notch (normalisee) 
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DFT du signal filtre: avec et sans la 
partie transitoire du signal filtre 
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Normalized Frequency ( *n red/sample) 
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Correlation de si avec si et s2 




TP #1 , Q2 



outl = conv (data , f liplr (signal_l) ) ; 
figure, plot ( (-62 : 9999) , outl) ; 




TP #1 , Q2 



outl = conv (data , f liplr (signal_2) ) ; 
figure, plot ( (-62 : 9999) , outl) ; 
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TP #1, Q2 - overlap-add 



% nb de points du signal de sortie: 1000+63-1 
% coordonnees absolues : -62:9999 
% coordonnees matlab: 1:10062 

outl = zeros (1,10000+63-1) ; %vecteur de 10062 zeros 
% signal decoupe en segments de 1000 
for i=0 : 9 

data_s = data(i*1000 +1: (i+l)*1000); 
out_s = conv (data_s , f liplr (signal_l) ) ; 
% data_s: i*1000 +1: (i+l)*1000; 
% f liplr (signal_l) : -63:-l 

% out_s : +63 - diff en coord, absolues et matlab 
outl (-63+1*1000+1+63: (i+1) *1000-l+63) = 

outl (-63+1*1000+1+63 : (i+1) *1000-l+63) + out_s ; 

end 

figure, plot ( (-62 : 9999) , outl) ; 
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TP #1 , Q3 

% filtre inverse 
rii_a = rif 1 ; 
rii_b = 1 ; 

% pour obtenir les valeurs des zeros de rifl 
abs (roots (rifl) ) ; angle (roots (rifl) ) ; 

% filtrage d'une colonne de 1' image 
img_Yr ( : , i) = filter (rii_b , rii_a , imgl ( : , i) ) ; 
% affichage de 1' image 
figure, imshow(img Yr, [0 255]) 
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Poles et zeros - filtre RIF1 
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Reponse en frequence du filtre RIF1 
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Poles et zeros - f iltre inverse 
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Reponse en frequence du f iltre inverse 





Image sortie 



Poles et zeros - f iltre RIF2 
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Filtre inverse instable 
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Une colonne de I' image 
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TP #1 , Q3b) 



% coefficients "a" du filtre inverse 
aa = rif 2 ; 

% equation aux diff . a resoudre 

% a(0)y(n) + a(l)y(n-l) + a(2)y(n-2) + a(3)y(n-3) + 
a(4)y (n-4)+a(5)y (n-5) = delta(n) 

% isole y(n-5) 

% y(n-5) = -(l/a(5)) (a(O)y(n) + a(l)y(n-l) + a(2)y(n-2) 
+ a(3)y(n-3) + a(4)y(n-4)) + delta (n) /a (5) ; 

% premier point non nul : y(-5); a(i) -> aa(i+l) 

yy = zeros (256,1) ; % n=0 -> yy(256); n=-255 -> yy(l) 

yy(-5+ 256) = l/aa(6); % solution pour n=0 
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TP #1 , Q3b) 



% equation aux diff . a resoudre 

% y(n-5) = -(l/a(5)) (a(O)y(n) + a(l)y(n-l) + a(2)y(n-2) 
+ a(3)y(n-3) + a(4)y(n-4)) + delta (n) /a (5) ; 

for n=-l : -1 : -249 

yy(n-5 + 256) = - (1/aa (6) ) * (aa (1) *yy (n+256) + 
aa(2)*yy(n-l+256) + aa (3) *yy (n-2+256) + aa (4) *yy (n-3+ 
256) + aa(5)*yy(n-4+256) ) ; 

end 

figure , stem(yy) ; 

[Li, Co] = size(img2); 
for i=l:Co 

img_Yr3 ( : , i) = conv (img2 (: ,i) ,yy) ; 

end 

figure, imshow(img Yr3,[0 255]) 
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Reponse a I'impulsion du f iltre stable 
256 premiers points avec retard de 256 
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Image filtree: filtre 
RIF de 256 pts 
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Image filtree: filtre 
RIF de 64 pts 
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Reponses en frequence des f litres 
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